Atomistic Coupling between Magnetization and Lattice Dynamics from First Principles 
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We formulate a new scheme to study the combined magnetization and lattice dynamics in magnets, so- 
called magnetoelastics. The coupling between magnetization and lattice are considered through an expansion 
of electron-phonon coupling, while the magnetization is coupled to electrons through sd-\ike interaction. We 
show the that the time-scale of the magnetization dynamics due to coupling between magnetic degrees of free- 
dom and electronic degrees of freedom can be transferred to lattice degrees of freedom and therefore can give 
rise to lattice dynamics on the same time-scale. This opens a new route to understand and treat ultrafast lattice 
dynamics induced by magnetization dynamics. We also show that all the parameters necessary to simulate this 
coupled lattice and magnetization dynamics can be obtained from first-principles. 

PACS numbers: 75.80.+q, 72.25.Rb, 71.70.Ej, 75.78.-n 



Magnetic phenomena at very fast time-scale is already a 
hot topic of both theoretical HI 01 as well as experimental re- 



search [|3|-|5|]. Magnetization dynamics at very fast scale is a 
primary requirement for efficient data storage since the writ- 
ing the data in magnetic way involves the process of switch- 
ing or magnetization reversal. Such memory devices such as 
magnetic random access memory (MRAM) is proposed to be 
a universal memory device fl. Similarly there are proposal as 
well as demonstrations of ferroelectric random access mem- 
ory devices (FeRAM) in which electric polarization due to 
mainly ionic displacements is considered as a variable for the 
construction of a bit. 

With the advances in the synthesis as well as theoretical 
understanding in multiferroic materials one goal of current re- 
search is aiming towards a possible memory device where ad- 
vantages of FeRAM and MRAM will be combined In 
the magnetic and ferroelectric case, bits are constructed from 
two states i.e, (M, — M) in magnetic case and (P, — P) in fer- 
roelectric case (M and P being the magnetization and electric 
polarization, respectively). In contrast, multiferroics materi- 
als offer multiple states such as (+M, +P), (+M, -P), (-M, 
+P), and (— M, — P) can be considered. One can, moreover, 
write data electrically and read magnetically, thereby reducing 
the energy consumption. 

It has been very difficult to realize dynamics in such sys- 
tems since polarization dynamics in many materials is much 
slower compared to the magnetization dynamics. The time 
evolution of electric polarization is related to lattice dynam- 
ics which can be very slow (about nano seconds) compared to 
the very fast magnetization dynamics (in pico or even femto 
second scale). Recently, however, Walker et. al have demon- 
strated experimentally lattice dynamics in femto-second scale 
in the multiferroic material TbMn03 |@], showing magnetic 
field induced ionic displacements. Here one should note that 
ferro-electricity in TbMnC>3 otherwise is not due to ionic dis- 
placement as in a proper ferroelectric like BaTi03. Cou- 
pling between lattice and magnetization has also been re- 
ported to have strong effects for iron based superconductors 
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SmFeAsOi-jFv and Bai ^K r FeAs 2 OH, LaFeAsOF 
1211 . and BaFe2As2 111 ll 1 1 311 ^ and in spin-ladder systems 



e.g. BaFe2Se3 01411 . It is important to notice that those experi- 
ments open a path-way to manipulate the polarization dynam- 
ics through the magnetization dynamics in such materials on 
a very short time-scale. 

We here develop a theoretical formalism which can be used 
to study ultrafast lattice dynamics in magnetic materials where 
magnetization is coupled to the ionic displacements. There 
are very few attempts in the literature where magnetization 
and lattice dynamics are studied as a coupled dynamical sys- 
tem. Ma et. al II 1511 had developed a model for simulat- 
ing molecular dynamics for magnetic material where a co- 
ordinate dependent exchange interaction between the spins 
determines the coupling between the spin and lattice degrees 
of freedom. The co-ordinate dependent exchange function is 
obtained by minimizing the expectation value of energy of a 
Stoner Hamiltonian with respect to magnitude of the magnetic 
moments. One should remember here that such minimization 
strongly neglects the spatial variation of charge density in the 
material and such approximation is valid only for very good 
metal where local c harg e neutrality is a natural situation. Re- 
cently Wang et. al 11611 also have developed a scheme based 
on molecular dynamics, where ionic dynamics was obtained 
from Newton's equation of motion while magnetization dy- 
namics was obtained from stochastic Landau-Lifshitz-Gilbert 
(LLG) formalism II171 1 1 811 . The authors had studied complex 



electric and magnetic susceptibilities of multiferroic BiFe03 
using such method. 

The LLG formalism is a restricted approach to study the 
magnetization dynamics in the sense that the coupling be- 
tween magnetic degrees of freedom with environmental ones 
are not considered in a proper way. The means of dissipation 
considered is only a single parameter: the Gilbert damping 
parameter. In recent developments [ 19 20l very fast dynam- 



ics have be demonstrated by coupling the magnetization with 
surroundings. Aiming at describing magnetism dependent lat- 
tice dynamics, we here show that one can obtain an effective 
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coupling between the magnetization and lattice dynamics, by 
inclusion of charge-lattice coupling. We derive coupled equa- 
tions for magnetization and lattice dynamics which can be 
implemented in existing first-principles codes, and all the pa- 
rameters needed for the simulation can be obtained from the 
first principles. Thus, the purpose of this letter is therefore 
two-fold: (i) design a scheme to perform combined magne- 
tization and lattice dynamics which goes beyonds methods 
based on Newtonian mechanics combined with classical LLG- 
equation, (ii) address ultrafast lattice dynamics observed re- 
cently in experiments. For the simulations in the microcanon- 
ical level (closed system), our method provides a direct way to 
study the energy transfer between different subsystems, such 
as magnetization, phonon and electrons. Also, our method 
can be used to study the recently observed multiglass behav- 
ior (both spin and dipole glass) in the some doped perovskites 
such as Mn doped SrTiC>3 12 111 . 

Our model for the magnetic interactions begins from the 
assumption that the magnetization M(r) interacts with the 
surrounding spin density s. s (r') through the potential v(r) = 
— f 7(r,r')cr -s^r^c/r'. We, thus, write the interaction Hamil- 
tonian Jf M = (1/2) / y^(r)v(r)v(r)dr = - f /(r,r')M(r) • 
s s (r')drdr', where ^/(r) = (i/^-(r) V / |( r )) r is a spinor such 
that we define M(r) = v/ t (r)en//(r)/2, whereas J(r,r') = 
J(r',r). 

The charge n(r) = So(r) = y/^ (r) y/(r) is subject to the po- 
tential 0(r) = / 0(r,Q(r'))<ir' due to electron-ion interac- 
tions, where Q(r) is the ionic displacement from its equilib- 
rium position. Here, we do not assign any specific nature of 
the displacement. However, in the perspective of multiferroic 
materials, for instance, one can think of Q(r) representing a 
polar mode or a mode which corresponds to rotation of oxy- 
gen octahedra. For small displacements, we employ the ex- 
pansion (r , Q (r ' ) ) w 0o (r ) + Q (r' ) • V r / O (r ) , where O (r ) = 
0(r,O), which gives the interaction between the charge and 
lattice vibrations = / Q(r') • V r /0o(r)«(r)<irc/r / . 

Given the general non-equilibrium conditions in the system, 
e.g. temporal fluctuations and currents, we define the action 



namics, is given by (using x = (r,f )) 
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(1) 



on the Keldysh contour J22-24]. Here, — 
— y§ f B ex t(r,f) ■M(r,t)drdt represents the Zeeman coupling 
to the external magnetic field B ex t, whereas S^wzwn = 
J§ /oM(r,f;T) • [<? T M(r,f;T) x d,M(r,t; x)]dxdt |M(r)| 2 Jr 
describes the Berry phase accumulated by the spin. The free 
lattice is represented by ,9\ att = § f{[/Q(r,f) • d t Q(r,t) — 
p 2 /2M ion ]<5(r-r') -Q(r,f) -V rr /(f) -Q(r',t)}drdr'dt, with 
the ionic mass M; on and the dyad V rr / = V r (V r /Vo), and 
where Vo is the ionic potential at equilibrium (vanishing 
displacements). 

Our model for the coupled magnetization and lattice dy- 
namics is obtained as a second order cumulant expansion of 
the partition function J°[M(r,f ),Q(r,f )] = tr T c e ir , tracing 
out the electronic degrees of freedom. Thus, the effective ac- 
tion S^mp, which describes the magnetization and lattice dy- 
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Q(x) • 3> m (x,x') -Q0O+M(jc) • -Q(x') 



+ Q(jc) ■ ® 0s (x,x') -M(x) +M(jc) • % s {x,x) -M(x )jdxdx, 

(2) 

where we have defined 



%q{x,x') = j E p {r,p)K pq {p,p';t,t)E q {p' ,r')dpdp' , 



(3) 



with the notation E (r,p) = V r o (p), E,(r,p) = -J(r,p), 
and K pq {p,p';t,t') = (-i)(Ts p (p,t)s q (p,t)), p,q = 0,s. 

The effective action contains phonon-phonon, spin-phonon, 
and spin-spin interactions, represented by the first, second and 
third, and fourth terms, respectively. The phonon-phonon in- 
teraction account for anharmonic effects in the material and 
the interactions is mediated by electronic density-density cor- 
relations (^oo)- The spin-phonon part is mediated by cor- 
relations between the spin and electron density (@ s o, ^0s)- 
For collinear magnetic background, this interaction reduces 
to the correlations between the electronic density and the spin- 
polarization acting solely in the z-direction. For chiral mag- 
netic structures, e.g. spin-orbit coupled systems, however, the 
interaction has finite components in all spatial directions. It 
is, hence, clear that the magnetic structure of the electronic 
background may have a significant influence on the mechan- 
ical subsystem. The electron mediated exchange interaction 
&n was discussed in 1201. 



It is worth mentioning that the magnetomechanical, or mag- 
netoelastic, coupling not only depends on the background 
electronic structure, but also to large extent on the coupling 
/ between the local spins and the electronic structure, on the 
one hand, and the potential landscape 0(r) and its gradient 
V0 experienced by the charge, on the other hand. It is, thus, 
easily understood that materials with slowly varying potential 
0, such that V0 « 0, are expected to have weak or negligible 
magnetomechanical coupling. This would be pertinent to ma- 
terials with no or uniform displacement distribution, e.g. fer- 
roelectric materials. This should be contrasted with magnetic 
materials in which a spatially fluctuating ionic potential leads 
to a non-negligible interaction parameter V0 which, in turn, 
implies a significant magnetomechanical coupling, pertaining 
to e.g. TbMn0 3 |9] and BaFe 2 As 2 fiUl . 



The equations of motion for the magnetization M and dis- 
placement Q are found by variation of the total action S^r = 
•S^iatt + ^mp + ^wznw + y<!xt with respect to fast fluctuations, 
see e.g. Ref. ll20ll for details. Requiring d t |M(x) | 2 =0, we ob- 
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tain 



M(x) =M(jc) x -yB ext (x) 



Ux,x') ■ Q(x') + &„{x,tI) • M(x') j tbi 

(4a) 



M i0 „Q(x) =- J (V^ -Q(r',*) + Q(r',t) ■ V rV J dt> 

+ \j (% S M ■ M(x') + <%oM • Q(*')) dx'. 

(4b) 

The equations for M and Q provide a general framework 
for coupled treatment of magnetization and lattice dynam- 
ics. It is legible that time evolution of both displacement and 
magnetization depend on surroundings in terms of space and 
time. Thus, Eq. emphasize that both Q and M, in gen- 
eral, depend non-locally on both the magnetization and ionic 
displacements for the entire structure. This is in sharp con- 
trast to previous exisiting methods which normally are based 
on molecular dynamics where Q and M are derived from ef- 
fective total energy in the adiabatic approximation, i.e. instan- 
taneous arrangements of magnetization and positions of ions. 



J 



In this context it is interesting to observe that the time evo- 
lution of a local mode B25I1 . is non-locally influenced by the 
magnetization at different points in space and time. It can, 
moreover, be seen that the ionic dynamics can be controlled 
by external fields, e.g B ext (r,f), something that was experi- 
mentally demonstrated in the Ref . J^] . 

The non-locality prescribed in the equations of motion 
for the magnetization and displacement, Eq. ©, provides 
a firm foundation for a coupled dynamics of the two vari- 
ables. For a deeper insight into the process that dominate 
the dynamics of the two variables, we here outline an ap- 
proximation scheme which is essentially equivalent to the 
Landau-Lifshitz-Gilbert paradigm of magnetization dynam- 
ics. Thus, considering the magnetization and lattice dynamics 
to be slow compared to the electronic processes, we expand 
in time M(r',f') ps M(r',f) - TM(r',f) + T 2 M(r',f)/2 and 
Q(r'/) ps Q(r,f ) - rQ(r',f) + T 2 Q(r',f )/2, where x = t-t'. 
We define the internal magnetic field B; nt (x) = — J @ ss {x,x!) • 
M(r' ,t)dx' / j - J & s q{x 1 x') ■ Q(r',f )dx'/y such that the to- 
tal magnetic field B(x) = B ex t(^) +Bi nt (x), and the internal 
force field F(x) = J(V rr , • Q(r',f) + Q(r',f) • V r , r )dr'/2 + 
J[% s (x,x') • M(r',f) + % (x,x') ■ Q(r' ,t)]dx' /2. Then, we 
can write the resulting equations of motion for M and Q as 



M(x) =M(jc) x (-yB(x) + |[G.„(x,rO-M(^0+U*y)^^ (5a) 

(5b) 

I 



M ion Q(x) =f(x)+- /[c^(jcy) -0(^,0 -0(^,0 +<^(*.o-^»o+w*y)-*(^oK 
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In Eq. §5^ we have introduced the intrinsic mag- 
netic (mechanical) damping and moment of inertia ten- 
sors G ssm (x,r r ) = - fT@' ssm (x,x')dt' and I ii(00) (x,r') 

j T 2 ^s(oo)(x ,x')dx/2, respectively, which appear also in any 
uncoupled description of the magnetization and lattice dy- 
namics. In addition to those intrinsic interactions, Eq. (0 
also contain corresponding magneto-mechanical damping and 
moment of inertia tensors G pq (x,r') = —J T2? pq (x,x')dx' and 
r') = / x 2 S l pq (x,x')dx/2, thus, showing that the mag- 
netic (mechanical) subsystem has a strong influence on the 
mechanical (magnetic) one. It, moreover, suggests that impor- 
tant time-scales for the two subsystems are governed by anal- 
ogous electronic processes, since the propagators K pq (x,x') 
all have similar dependencies on the electronic structure. We 
can, therefore, infer that the dynamical processes of the mag- 
netic and lattice degrees of freedom occur on the same time- 
scales, from which we can also understand that dynamical 
variations in the lattice (magnetic) subsystem transfer to the 
magnetic (lattice) subsystem without significant delay. We 
expect that this conclusion will be of particular importance 



for understanding the magneto-mechanical dynamics that has 
been observed in e.g. multiferroic and superconducting mate- 
rials. 

Next we describe how the parameters in Eq. (0 can 
be calculated from a first principles point of view. The 
conditions defined by the DFT system, the time local- 
ity of & pq enables the Fourier transform ^ r (r,r';e) = 
/ 2>(x,x')e ,£T dx' . Thus, we can write the damping and iner- 
tial tensors in terms lim £ ^o /<9 £ £^ (r,r';e) = / x@l q {x,x!)dt' 
and lim £ ^o<5 £ / X 2 £t pq (x,x')dt', where 

%Jry-e) =4sp / Z p (r,p)*PlmG r j((D)*«hnG r p A(D>) 



„ , , ,. f(<a')-f(co) dcodco' , 



£ - CO + CO' 



i8 2k 2k 



(6) 



where G' rr ,(co) is the retarded Green function for the back- 
ground electronic structure represented as a 2 x 2 matrix in 
spin space over which the trace sp is taken. Here also p,q = 
0,s using er° — <7° (2 x 2 identity) and a" = er. We notice 
that this expression presents a generalized form of exchange 



4 



interaction in both magnetic and lattice subsystems as well as 
between the two of them. Using Kramer's-Kronig relations in 
the limit e — > we can write this generalized exchange inter- 
action as 



1 



(r/)=-spIm//((D)E p (r,pKGp p ( 



x(r e 'G' p ,(co')Z,(p',r')dcodpdp' (7) 



'pq\ 



in accordance with Ref. [20] for the purely magnetic case. 
The generalized damping tensor can be written as 

G M (r,r') =lsp J /'(©)S p (r,p)^ImG^ p (fi)) 

x <TnmG> ppl {co)Z q {p' y)dcodpdp' . (8) 

In this form, it is easy to see that the main influence from 
the electronic structure on both the magnetic and mechani- 
cal damping, as well as the magneto-mechanical damping is 
generated by the properties at the Fermi surface. In the same 
fashion we can write the generalized moment of inertia tensor 
according to 

I w (r,r') =sp J /(ffl)S p (r,p)pImG p , p (G)) 

x^ R eG pp ,(®) + ^ImG pp ,(co) 

x o-^ReG^O))] E q (p',r')^dpdp', (9) 

showing that the moment of inertia depends on the occupied 
portion of the electronic structure. The derivatives of the GFs 
illustrate that the moment of inertia is very sensitive to sharp 
variations in the electronic structure. We, therefore, expect 
that the moment of inertia may have significant effect in nar- 
row band magnetic materials, e.g. strongly correlated electron 
systems. Moreover, in a similar fashion as the magnetic order 
may strongly influence the generalized exchange it may also 
have a strong influence on the generalized moment of inertia. 
The coupling between the magnetic and mechanical subsys- 
tems indicates that the magnetic order also has a significant 
influence on the mechanical systems. In an analogous way 
we, furthermore, expect the lattice ordering to have a signifi- 
ant effect both on the mechanical as well as on the magnetic 
dynamics. 

Despite the apparent similarities of Eqs. © — © with the 
corresponding equations provided in 112011 . the purpose of pre- 
senting them here is to show that not only the magnetic inter- 
actions, particularly the exchange and damping, can be given 
in a form of Kubo formula, but that this also holds for the me- 
chanical and magneto-mechanical interactions. This is very 
appealing from the perspective of the fluctuation-dissipation 
theorem as it demonstrates that the coupling between the mag- 
netic (mechanical) fluctuations and mechanical (magnetic) 
dissipation, can be treated on the same footing as the pure 
magnetic and mechanical correspondences. This very clearly 
illustrates the important role of the charge density mediated 
interactions, both for the magnetic and mechanical systems, 
as well as the interactions between them. 



The above two coupled equations, Eq. (01 or Eq. (0, pro- 
vide a tool for simulating the dynamics of a closed system 
composed of electron, magnetization, and phonon subsystems 
and various time-resolved quantities such as spin-spin or elec- 
tric dipole-dipole autocorrelation functions can obtained eas- 
ily. Also, one very interesting feature of our formulation is 
that, this in principle should work for the cases where spa- 
tial distribution of charge density is quite non-uniform. We 
have explicitly considered the charge densities in terms of the 
correlation function K in our formulation. Further extension 
of this work would be to develop a scheme for controlling the 
temperature of the system during the course simulation, which 
is beyond the scope of present study. 
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